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Abstract 

We study spectra, localization properties and local chirality of eigen- 
vectors of the lattice Dirac operator. We analyze ensembles of quenched 
SU(3) configurations on both sides of the QCD phase transition. Our 
Dirac operator is a systematic expansion in path length of a solution of 
the Ginsparg- Wilson equation. Analyzing the finite volume behavior 
of our observables and their scaling with the gauge coupling we come 
up with a consistent picture of topological excitations on both sides of 
the QCD phase transition. Our results support the interpretation of 
the dominant gauge field excitations seen by the Dirac operator as a 
fluid of instanton-like objects. 
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1 Introduction 



The QCD phase transition tightly knits together two of the most intriguing 
features of QCD: When the temperature drops below the critical value Tc, 
chiral symmetry is broken and quarks become confined. State of the art 
lattice results indicate that both these transitions happen at the same 
(see e.g. |l| for a recent comprehensive review of finite temperature lattice 
results). Lattice studies which compare the changing behavior of observ- 
ables on both sides of Tc might help to better understand the underlying 
excitations leading to confinement and chiral symmetry breaking. 

While different ideas for excitations responsible for confinement are still 
heavily debated, the understanding of chiral symmetry breaking seems to 
be in better shape and during the last 20 years a favorite candidate for 
the relevant degrees of freedom has emerged. It is now widely believed 
that models of instanton-like excitations provide a good description of chiral 
symmetry breaking (for extended reviews see e.g. These models have 

a wealth of observable signatures and here we set out to test these predictions 
in an ab initio lattice calculation and we compare the observables on the two 
sides of the QCD phase transition. 

The basic idea of chiral symmetry breaking from instantons is based on 
two arguments: Firstly, the Banks-Casher relation |Q] (see Eq. ( |3.lD below) 
relates the chiral condensate ("0^) to the density p(A) of eigenvalues A of 
the Dirac operator near to zero. Secondly, the near zero modes which build 
up the density of eigenvalues p(0) at small A are believed to come from 
interacting instantons and anti- instantons: As long as an instanton is well 
separated from an anti-instanton, the Dirac operator produces only zero 
modes which do not contribute to p{Q) in the Banks-Casher formula. As soon 
as the instanton and an anti-instanton start to overlap, the Dirac operator 
can no longer distinguish them as isolated topological objects and instead of 
two zero modes produces a pair of complex conjugate eigenvalues. For small 
overlap, these two eigenvalues are still very close to the origin but move 
up (down) the imaginary axis as the instanton and anti-instanton approach 
each other. 

While earlier studies of instantons on the lattice mainly were based on 
direct filtering of the gauge fields such as cooling the recent progress 
for chirally symmetric lattice fermions based on the Ginsparg- Wilson re- 
lation Q now allows to directly use the spectrum and eigenvectors of a 
chirally symmetric Dirac operator to test the instanton picture (see 0"[2C] 
for studies of spectrum and eigenvectors of different lattice Dirac opera- 
tors). When using the eigenvectors of the Dirac operator as a probing tool 
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for instanton-like excitations, the basic assumption is that the locahzation 
of the eigenvector traces the locahzation of the underlying instanton config- 
uration, a property which is known for eigenvectors with eigenvalue of the 
continuum Dirac operator in an instanton or multi-instanton background. 
The scenario for chiral symmetry breaking sketched in the last paragraph 
leads to many characteristic effects that should be visible in the spectrum 
and eigenvectors of the Dirac operator. In particular the spectral density, 
the localization of eigenvectors and their local chirality give characteristic 
signatures which can be tested in an ab initio calculation on the lattice and 
which might support the instanton picture or rule it out. 

In this article we report on our lattice studies of eigenvectors and eigen- 
values of the Dirac operator. We use ensembles of SU(3) gauge fields in the 
quenched approximation. Our runs were done on three different lattice sizes 
in order to study the volume dependence of observables. The gauge cou- 
plings were chosen such that we have ensembles on both sides of the QCD 
phase transition. This allows to directly compare the changing behavior 
of observables as one changes from the confined, chirally broken phase to 
the deconfining, chirally symmetric phase. Such a direct comparison reveals 
which degrees of freedom show the most drastic changes and thus are rel- 
evant excitations for chiral symmetry breaking. Our Dirac operator is an 
approximate solution of the Ginsparg- Wilson equation based on a system- 
atic expansion [pl|]- |]2^ and we will refer to it as the chirally improved Dirac 
operator. It is considerably cheaper to implement than e.g. the overlap op- 
erator [p^ j while producing very similar results for the relevant eigenvectors 
and eigenvalues in the physical branch ||2^ . This property makes the chirally 
improved Dirac operator ideal for the present study since it allows for very 
large statistics which would not have been possible for the overlap operator. 

We simulate finite temperature QCD by working on a lattice with Lt < 
Lg, using periodic time boundary conditions for the gauge fields, and an- 
tiperiodic for the fermions. The temperature is given by T = l/Lf. 

Let us finally remark on the language convention used in this paper: 
The short and periodic time-direction Lt introduces an additional length 
scale which is felt by extended excitations of the gauge field. The solution 
of the SU(3) equations of motion at finite temperature is usually referred 
to as caloron. Calorons exhibit a string-like localization in time direction. 
However, the transition between calorons and instantons is continuous, since 
very small instantons hardly feel the finite time extent. Thus we here use 
the term topological excitations and speak of instantons or calorons only 
when we want to stress particular differences. 

The article is organized as follows: In the next section we discuss the 
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setting of our computations which encompasses a brief presentation of the 
chirally improved Dirac operator (Section 2.1), a discussion of the gauge 
ensembles (Section 2.2) and the techniques used for the diagonahzation of 
the Dirac operator (Section 2.3). In Section 3 we collect our results for the 
spectrum of the Dirac operator. These results are presented in two parts, 
Section 3.1 where we discuss the density distribution and Section 3.2 where 
we analyze the fermionic definition of the topological charge via the index 
theorem. In Section 4 we present our results for the eigenvectors. In Section 
4.1 we give a qualitative discussion based on some examples followed by a 
discussion of the localization properties of the eigenvectors in Section 4.2. In 
Section 4.3 we present our results for the local chirality of the eigenvectors. 
The paper closes with a summary in Section 5. In two appendices we give 
the details of our Dirac operator (Appendix A) and discuss the fermion zero 
modes in a classical caloron solution (Appendix B). 



2 Preliminaries 



2.1 The chirally improved Dirac operator 

In the current analysis we use a chirally improved Dirac operator D which 
is an approximation to a solution of the Ginsparg- Wilson equation D 
is a systematic expansion of a Ginsparg- Wilson Dirac operator in terms of 
paths on the lattice. A detailed description of this approach was presented 



in ||2l|, 23 1 and here we only sketch a few basic steps of the construction. 

For about three years it has been understood, that chiral symmetry can 
be implemented through a Dirac operator D which obeys the Ginsparg- 
Wilson equation @ (we set the lattice spacing to 1) 

75 I? + Dj5 = Dj^D. (2.1) 



Currently two types of solutions are known, the so called overlap operator 
[p4| and perfect actions ^] . The overlap operator gives a simple expres- 
sion for a Dirac operator obeying Eq. ( |2.1| ), but the involved square root of a 
large matrix, makes a numerical evaluation quite costly. In two dimensions 
the fixed point Dirac operator has been thoroughly tested in the Schwinger 
model In four dimensions the perfect Dirac operator has entered the 
test phase |2^, but so far no explicit parametrization has been published. 

A third approach to the Ginsparg- Wilson equation is the above men- 
tioned systematic expansion of a solution of Eq. ( |2.1| ). The first step in 
the construction of our chirally improved Dirac operator is to write down 
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the most general Dirac operator on the lattice. This is done by allowing 
more general lattice discretizations of the derivative. The standard deriva- 
tive term makes use of nearest neighbors only but certainly one can include 
also more remote points on the lattice such as next to nearest neighbors or 
diagonal terms etc. Each such term is characterized by the product of link 
variables which form the gauge transporter connecting the two points used 
in the derivative. The corresponding set of links can be viewed as a path 
on the lattice. The most general derivative on the lattice will then include 
all possible paths, each of them with some complex coefficient. In order to 
remove the doublers, in addition to the derivative terms coming with the 
Dirac matrices 7^ we also have to include terms proportional to the unit 
matrix in Dirac space. To obtain the most general expression we include all 
16 elements Fq, of the Clifford algebra, i.e. we also add tensor, pseudovector 
and pseudoscalar terms. To summarize, our Dirac operator is a sum over all 
Fq,, each of them multiplied with all possible paths on the lattice and each 
path comes with its own coefficient. 

The next step is to apply the symmetry transformations: translations, 
rotations, charge conjugation, parity and 75-hermiticity which is defined by 



are used to prove CPT in the continuum. Once these symmetries are imple- 
mented the coefficients of the paths in the Dirac operator are restricted. One 
finds, that groups of paths which are related by symmetry transformations 
have to come with the same coefficient, up to possible signs. We denote the 
most general Dirac operator which obeys the symmetries as: 



^75 



(2.2) 




can be seen to correspond to the properties of D which 
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75 s{ll) sih) s{k) s{k) pi < li,lj,lk,ln > 



ii=±li'2=±2 ij,k,n=l 
i3=±3,i4 = ±4 



(2.3) 



By e we denote the totally anti-symmetric tensors with 2,3 and 4 indices. 
We use the notation < li,l2 ■■■In > to denote a path of length n and the 
li G {— 4, — 3, — 2, — 1, 1, 2, 3, 4} simply denote the directions of the subse- 
quent links which build up the path. s(/j) is an abbreviation for sign[li). 
With the particular choice for the generators of the Clifford algebra used in 
Eq. ( |2.3D (no additional factors of i), the coefficients Si,Vi,ti,ai,pi are real. 
The expansion parameter for the Dirac operator in Eq. ( p.3| ) is the length 
of the path since the coefficients in front of the paths decrease in size as 
the length of the corresponding path increases. A general argument for this 
behavior can be given and it has been confirmed numerically for the solution 
presented in ||2^. We remark, that an equivalent form of D presented in [^] 
is the basis for a parametrization of the perfect Dirac operator. 

The final step in our construction is to insert the general expression for D 
into the Ginsparg- Wilson equation. On the left hand side of the Ginsparg- 
Wilson equation ( |2.lD some of the terms acquire minus signs, depending 
on the commutator of the corresponding Ta with 75. On the right hand 
side an actual multiplication of all the terms in D has to be performed. 
However, using the above notion of a path, the multiplication on the right 
hand side can be formulated in a algebraic way and then can be evaluated 
using computer algebra. Once all multiplications are performed one can 
compare the left and right hand sides of the Ginsparg- Wilson equation. It 
is important to note that for an arbitrary gauge field different paths, which 
correspond to different gauge transporters, are linearly independent and can 
be viewed as elements of a basis. Thus for the two sides of Eq. ( p.l| ) to be 
equal, the coefficients in front of the same basis elements on the two sides 
have to agree. When comparing the terms on the two sides, the result is a set 
of coupled quadratic equations for the expansion coefficients Si,Vi,ti,ai,pi. 
This set of equations is equivalent to the Ginsparg- Wilson equation. After a 
suitable truncation of ( |2.3| ) to finitely many terms the system can be solved 
and the result is an approximation to a solution of Eq. (^]^) . In addition it 
is possible to include a dependence on the inverse gauge coupling /3i through 
an additional constraint for the coefficients. This step allows to work with 
less terms in the parametrization. This procedure is similar to the tuning of 
the mass-like shift which is used to optimize the localization of the overlap 



operator |3C]. An explicit list of the terms used in our parametrization of 
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the Dirac operator and the values of the coefficients Si,Vi,ti,ai,pi are given 
in Appendix A. 

After a test of the 2-d chirally improved Dirac operator in the Schwinger 
model with dynamical quarks in ||2^ the construction was outlined for four 
dimensions in [^]. A test of a Dirac operator based on this approximation 



was presented in [23| and it was found that the approximation is particularly 
good in the physical part of the spectrum. It was found that near the origin 
the deviation of the eigenvalues from the Ginsparg- Wilson circle is very 
small. Since in the present study we are interested in the spectrum near the 
origin, this makes the chirally improved Dirac operator very well suited for 
the physical questions analyzed here. It is numerically much less demanding 
than the overlap operator and has a well ordered spectrum near the origin 
such that spectral densities and the formation of a gap can be studied. 

Before continuing with the discussion of the gauge action in the next 
subsection, let us discuss a general property of the eigenvectors of a Dirac 
operator D which obeys the 75-hermiticity of Eq. (p.2|). Let if; be an eigen- 



vector of D, i.e. Dip = X^p. Then a few lines of algebra show |31] 



"0^ 75 V' = ! unless A is real . (2-4) 

This equation is the basis for the identification of the zero modes. For the 
continuum Dirac operator or also for a D which is an exact solution of the 
Ginsparg- Wilson equation only zero modes have a non-vanishing matrix el- 
ement with 75. Here we are working with an approximate solution of the 
Ginsparg- Wilson equation and we do not have exact zero modes. However, 
eigenvectors where the corresponding eigenvalue has a non-vanishing imag- 
inary part are excluded as candidates for zero modes, since their matrix 
elements with 75 vanish exactly. In our approximation the zero modes show 
up as eigenvectors with a small real eigenvalue and a value of ip'^^/^ip ~ ±0.87. 
If one adds more terms in our approximation then the zero modes would be 
closer to and the matrix element with 75 closer to ±1. 

2.2 The gauge ensembles 

The ensembles of gauge configurations for our study were generated with the 



Liischer-Weisz action |^, |3^. To be specific, we use the improved gauge 



action as presented in [^]. Explicitly, the gauge field action reads 
S[U] = /3i 5] ^ Re Tr [1 - Upi] + /32 E ^ Tr [1 - Urt] 

pi rt 
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/33E^ReTr[l-C/g 

pg 



(2.5) 



where the first sum is over all plaquettes, the second sum over all 2 x 1 
rectangles and the last sum runs over all parallelograms. /3i is the principal 
parameter while /52 and /Js can be computed from (3i using tadpole improved 



perturbation theory |Q, giving |l33|| 



with 



P2 



no 



Pi 



20 ul 



[1 + 0.4805 q] 



Pi 

-i 0.03325 a 



(^Re Tr{Upi) 



1/4 



a 



In (iRe TriUpi) 



(2.6) 



(2.7) 



^3 ' 3.06839 
The couplings P2 , P3 are determined self-consistently from uq and a for a 
given Pi . In Table || we list the values of the Pi used for our ensembles and 
our results for the expectation value {Upi)/3 of the plaquette. 



Pi 


8.10 


8.20 


8.30 


8.45 




+0.6294 


+0.6370 


+0.6432 


+0.6521 


P2 


-0.5474 


-0.5499 


-0.5533 


-0.5577 


P3 


-0.0512 


-0.0502 


-0.0495 


-0.0483 



Table 1: Parameters for our ensembles of gauge field configurations. We list 
the values of the Pi and the expectation value {Upi)/3 of the plaquette. 

In |2^ it was noted that when comparing the expectation value {Upi)/3 
of the plaquette for the Liischer-Weisz action with the corresponding expec- 
tation value for the Wilson action one finds that the former gives rise to 
a value of {Upi)/3 considerably closer to 1 (such a comparison makes sense 
only after matching a physical scale such as e.g. the string tension). Thus 
the Liischer-Weisz action tends to suppress ultraviolet fluctuations and typ- 
ically one obtains better results for approximate Ginsparg- Wilson fermions 
[ p3[ as well as for the overlap Dirac operator p5| ]. 

For all runs presented here the time extent of the lattice was held fixed 
at Lf = 6. Our values of Pi were chosen such that for j3i = 8.10 the system 
is in the confining phase, while the other three values {Pi = 8.20, 8.30, 8.45) 
are in the deconfined phase. Based on a scaling analysis of the spectral gap 
[ p6[ , we estimate that the phase transition occurs at f3i ~ 8.18. 
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0.1 






0.0 


/■ 
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-0.1 0.0 0.1 



Figure 1: The Polyakov loop in the complex plane. Each plot shows 800 
measurements for lattice size 6 x 16^. 



In order to set the scale we computed the Sommer parameter []37| vq for 
Pi = 8.10,8.30 and /3i = 8.45 on 16^ lattices |3§]. We give our resuhs for 
ro/a together with the lattice spacing a (assuming vq = 0.5 fm) in Table |. 

In Fig. m we show scatter plots of the Polyakov loop [^9| (L) in the 
complex plane. For each of our samples we display 800 measurements of 
(L) on a 6 X 16^ lattice. For f3i = 8.10 the measurements of the Polyakov 
loop are concentrated at the origin indicating that the system is in the 
confined phase. The other three plots show the emerging of three disjoint 
clusters of expectation values characteristic for the deconfining phase. The 
invariance of the pattern under rotations by 27r/3 is due to the Z3 symmetry 
of the gauge action. The action for fermions does not display this symmetry 
and thus fermionic observables behave quite differently in the sector with 
real Polyakov loop when compared to the two complex sectors [40|. For 
some of the observables below we will divide our ensembles of gauge field 
configurations into two sectors according to the phase of the Polyakov loop. 
We will refer to the configurations where the Polyakov loop (L) has a phase 
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~ as the real sector, while configurations with phase ~ ib27r/3 wiU be 
referred to as the complex sector. 

Finahy we comment on the statistics for our runs. For two values of (3i, 
one in the confining phase = 8.10), the other one in the deconfined phase 

= 8.45), we performed runs at three different values of the spatial extent 
Lg of our lattice {Ls = 12, 16, 20) in order to study the finite volume behavior 
of the system. The temporal extent Lt was kept fixed at Lt = 6. To study 
the /3i-dependence in the deconfined phase we performed two more runs on 
the 6 X 16^ lattice at /?i = 8.20 and (3i = 8.30. In order to compensate for 
the larger fluctuations on the smaller volumes we analyzed larger samples 
for these cases. The updates were done with a mix of Metropolis and over- 
relaxation steps. We also include a random Z3 rotation in order to update 
this symmetry of the action. Our configurations are decorelated by roughly 
two integrated autocorelations times of the Polyakov loop. Table || gives an 
overview over the statistics for our ensembles. 





Pi = 8.10 Pi = 8.20 


Pi = 8.30 


Pi = 8.45 


ro/a 


3.94(16) 


4.67(13) 


5.00(5) 


a 


0.127(5) fm 


0.107(3) fm 


0.100(1) fm 


6 X 12^ 


1200 




1200 


6 X 16^ 


800 800 


800 


800 


6 X 20^ 


400 




400 



Table 2: Sommer parameter tq, lattice spacing a and statistics for our gauge 
field ensembles. 



2.3 Technical remarks 

For the computation of the eigenvalues and eigenvectors we used the implic- 
itly restarted Arnoldi method [41|. For each of our gauge field configurations 



we always computed 50 eigenvalues. The search criterion for the eigenvalues 
was their modulus, i.e. we computed eigenvalues in concentric circles around 
the origin until 50 were found. 

When comparing these eigenvalues for different system sizes, one has to 
take into account the different density of eigenvalues for different volumes. 
For larger volumes the Dirac operator has a higher density of eigenvalues 
than on smaller volumes such that when keeping the number of evaluated 
eigenvalues fixed at 50, the last eigenvalues found have larger imaginary 
parts for smaller lattices when compared to the results for larger systems. 
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To give an example, in the chirally broken phase {Pi = 8.10) the largest 
eigenvalues we found among the 50 computed eigenvalues on the smallest 
lattice (6 x 12^) typically had imaginary parts of size ~ 0.3 while for the 
largest lattice (6 x 20^ ) the largest eigenvalues reached only imaginary parts 
of ~ 0.14 (compare the density plots in the next section). Some of our 
observables are sensitive to these differences due to the technical cut-off of 
50 computed eigenvalues and for them we cut off our data at a common 
physical value for all volumes, i.e. at |Im A| = 0.14. 



3 Results for the spectrum 

3.1 Density distribution of the eigenvalues 

Along with the transition to the deconfined phase, high temperature QCD 
also restores chiral symmetry, leading to a vanishing chiral condensate (ipip). 
Since the chiral condensate is related to the density^ p(A) of eigenvalues A 
at A = via the Banks-Casher formula Q 

(V^V) = -7rp{0)V-\ (3.1) 

one expects quite a different behavior of the Dirac spectrum on the two sides 
of the transition. In the deconfined, chirally symmetric phase the density 
of eigenvalues has to vanish at the origin, up to some zero modes due to 
isolated topological excitations in the gauge fields. In the confined, chirally 
broken phase a non- vanishing density of eigenvalues extends all the way to 
A = 0. 

In Fig. § we show three characteristic spectra, each consisting of the 50 
eigenvalues closest to the origin on 6 x 20'^ lattices. For the confining, chirally 
broken phase (Pi = 8.10, left hand side of the plot) we find a non- vanishing 
density at the origin. When going to the other side of the transition (Pi = 
8.45, the two plots on the right hand side of Fig. Q) one finds that the 
spectrum has developed a gap at the origin. There could be some isolated 
zero modes as e.g. in the central plot but they do not contribute to the 



density p{X) in the Banks-Casher formula (3J.). As already discussed above, 
the spectrum behaves quite differently for gauge configurations in the real 
and complex sectors of the Polyakov loop. For the example with the real 
Polyakov loop (central plot) we find a gap roughly twice as big as for the 
case with a complex Polyakov loop shown in the plot at the right hand side. 



^Our density p(A) is normalized such that p(A) AA gives the number of eigenvalues in 
the interval AA. 
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Pi = 8.10 

<L> = 



0.40 



0.20 



E 0.00 



-0.20 



-0.40 ^ ^ ^ 

0.00 0.20 

Re A, 



Pi = 8.45 

<L> real <L> complex 



0.40 



0.20 



E 0.00 



-0.20 



-0.40 ^ ^ ^ ^ — ^ 

0.0 0.2 0.0 0.2 



Re A, 



Re A, 



Figure 2: Examples for spectra of the Dirac operator. We show the 50 
eigenvalues closest to the origin, represented by circles. All spectra are for 
6 X 20^ lattices. On the left hand side we plot the eigenvalues A in the 
complex plane for /?i = 8.10 (chirally broken phase). The two plots on the 
right hand side show spectra at f3i = 8.45 (chirally symmetric phase) for 
a gauge configuration with the Polyakov loop in the real sector as well as 
for a configuration with complex Polyakov loop. The full curve in the plots 
indicates the Ginsparg- Wilson circle. 



There has been some debate whether in the complex sector the gap forms at 
all [ 40 1 , but we find that for values of Pi sufficiently above the critical value 
the existence of the gap can be clearly established. A detailed analysis of 
the volume dependence and the scaling of the spectral gap will be presented 
elsewhere Q. 

By averaging over all configurations in each sample (compare Table ^) 
one can compute the density p{X) of eigenvalues. For the continuum Dirac 
operator the spectrum is restricted to the imaginary axis and the definition 
of the spectral density is straightforward. For our lattice Dirac operator the 
spectrum is located on the Ginsparg- Wilson circle, or to be more precise in 
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0.1 0.2 

Im A, 



0.3 



Figure 3: Histogram plots for the spectral density p{X) normalized with the 
inverse volume as a function of Im A. 



the vicinity of this circle since here we are dealing only with an approximate 
Ginsparg- Wilson operator. In order to define the spectral density we bin the 
eigenvalues with respect to their imaginary part and count the number of 
eigenvalues in each bin. Eigenvalues with vanishing imaginary part, i.e. the 
zero modes were left out. After dividing the count in each bin by the total 
number of eigenvalues we obtain the histograms for p{X) used in Fig. |3| below. 
Since we are only interested in the density p{X) in a region of A where the 
real parts of the eigenvalues on the Ginsparg- Wilson circle are small, binning 
with respect to the imaginary parts of A gives a good approximation of the 
continuum definition of the spectral density. 

In Fig. |3| we show (see also ||9|, |lO|, 11 , |l2| ) the density of eigenvalues A 
normalized by the inverse volume as a function of Im A (we display only the 
positive half of the curve, i.e. Im A > 0). The top plot gives our results for 
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the confined, chirally broken phase. It is obvious that after normalization 
by the inverse volume, the data for different lattice sizes essentially fall on a 
universal curve. Only the edge near the origin is more rounded for smaller 
system size in agreement with universal random matrix theory predictions 
[42, 43]. For the same reason also the curve for the largest system (6 x 20^) 



still shows a drop near the origin. Up to this finite size effect, the density 
remains non-zero down to A = 0, thus building up the chiral condensate 
according to the Banks-Casher formula ( p.ip . The picture is quite different 
in the deconfined, chirally symmetric phase. Again we find that the data for 
different volumes fall on a universal curve, but now the density drops to zero 
already at nonvanishing values of Im A. For one of our samples (6 x 16^) 
in the real sector (bottom left plot in Fig. ^) we find a very small signal in 
the vicinity of the origin. We believe that this is a statistically insignificant 



fluctuation which does not show up for the larger systems. In [11| a similar 
event on a 4 x 16^ lattice with Wilson gauge action at /? = 5.75 and overlap 
fermions was observed and interpreted as a possible sign for a non-vanishing 
fermion condensate in the deconfining phase of quenched QCD. 

As already discussed above, the spectral gap is considerably smaller for 
the sector of gauge configurations with complex Polyakov loop but still is 
clearly visible (compare also 



3.2 Fermionic definition of the topological charge 

In this section we now analyze the zero modes, which for our Dirac operator 
correspond to eigenvalues with vanishing imaginary part. These modes come 
from isolated topological excitations of the gauge field. In the continuum the 
eigenvectors corresponding to the zero modes can be chosen as eigenmodes 
of 75 with eigenvalues ±1. The numbers n+,n_ of these left-, respectively 
right-handed modes are related to the topological charge Q via the Atiyah- 
Singer index theorem |44| 



Q = n+ - n_ . (3.2) 

A proof of the theorem for the case of a Ginsparg- Wilson lattice Dirac 
operator can be found in [^6|. Here we work only with an approximate so- 
lution of the Ginsparg- Wilson equation and we must give a more detailed 
description of our fermionic definition of the topological charge correspond- 
ing to Eq. (113) • 

There are two main effects of working with only an approximation of a 
solution of Eq. ( |2.1| ): Firstly, the topological modes do not appear at exactly 
zero, but typically have a small, nonvanishing real part which increases as 
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Figure 4: Distribution of the topological charge. The filled histograms show 
our numerical results and the dotted curve represents the simple instanton 
anti-instanton model discussed in the text. 
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the corresponding topological excitation of the gauge field becomes more 
like a defect, i.e. shrinks in size to about a lattice spacing (compare [|25[| ). 
Secondly, the corresponding eigenvectors ip are not exact eigenstates of 75 
but the absolute value of V'WsV' is smaller than 1, typically between 0.8 and 
0.9 for our D. Thus we define n_|_ (n_) to be the number of eigenvectors with 
positive (negative) values of ip^j^ip (as noted above ip^j^ip vanishes for all 
eigenvalues with nonvanishing imaginary part (see Eq. (2.4)). In addition we 
take into account only real eigenvalues smaller than 0.15 in order to discard 
defects. We should remark, that when moving this cutoff to either 0.10 or 
0.20 our results as e.g. shown in Fig. ^ below do not move beyond the error 
bars. This implies that for our D defects give only a negligible contribution 
to the fermionic definition of the topological charge. 

In Fig. 1^ we show our results for the distribution of the topological 
charge. The filled histograms show our numerical data and the dotted curve 
is the distribution from a simple model of non-interacting, dilute instantons 
and anti-instantons. 

The basic assumption is that in a cell of size 5V the probability for finding 
an instanton or an anti-instanton is given by pj 6V/2 while the probability for 
finding an empty cell is 1—pi 5V. Here pi is the density of the instantons and 
anti-instantons which are assumed to be distributed independently. After 
collecting the combinatorial factors one finds that the probability for having 
n-(_ instantons and n_ anti-instantons is given by a multinomial distribution. 
This is easily converted to a distribution P{Q) for the topological charge Q 
as given in Eq. ( |3.2D by summing over n_|_ and n_ with the constraint 6{Q — 
n+ + n_). When going over to the Fourier transform it is straightforward 
to perform the limit 5V which then gives the integrand in Eq. ( |3.3| ) 
below. Inverting the Fourier transform gives the final result 



P(Q) = — [ ^ dk e'^^ e-P'^^^'^"^^^ = e-P'^ Iq{piV) ■ (3.3) 
27r Jo 

Here Iq denotes the modified Bessel function of order Q and V is the volume 
of the system. The expectation value (Q^) is then given by 



2e 



-piV 



h{piV) + ^hipiV) + 9/3(^7^)... = piV . (3.4) 



Here we used a standard summation formula for the modified Bessel func- 



tions. Both results (3^) and ( p^ ) depend only on the dimensionless prod- 
uct pi V . Thus as soon as one computes (Q^) the distribution P{Q) is 
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Figure 5: Volume dependence of (Q^)- The numerical data are represented 
by diamonds and we connect the symbols to guide the eye. Please note the 
different scale of the vertical axis in the two plots. 



entirely fixed. We remark that in this model the topological susceptibility 
X = {Q"^) /V is simply the density pi of instantons and anti-instantons, i.e. 
X = PI- 

A discussion of more sophisticated dilute instanton models, which take 
interactions into account, can be found in e.g. [^]. 

When comparing our numerical results to the theoretical prediction (|3.3| ) 
in Fig. ^ we proceeded as follows: (Q^) was first computed for each combi- 
nation of f3i and volume (see our results below) and then used to evaluate 
the theoretical distribution P{Q). Thus the curves shown in Fig. ^ con- 
tain no fit. We find that on both sides of the transition and for all system 
sizes which we analyzed, the numerical data for P{Q) is quite accurately 
described by the simple instanton anti-instanton model. 

Equation ( [j.4D for {Q"^) also makes a simple prediction for the volume 
dependence of (Q^)- In Fig- § show our results for (Q^) as a function 
of the volume. For both phases the volume dependence is roughly linear. 
For the confined, chirally broken phase {(3i = 8.10, left hand side plot) one 
finds a considerably larger value of the topological susceptibility, i.e. the 
slope in the plot of (Q^) versus V (note the different scales for the two 
plots). This corresponds to a larger density pi in our model. It seems that 
in the chirally broken phase the assumption of diluteness for the instantons 
and anti-instantons becomes violated and the topological excitations can no 
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longer be viewed as non-interacting. This may be the reason for the sUght 
deviation from a strictly linear behavior of the f3i = 8.10 data. For a more 
accurate modeling of the Pi = 8.10 data the simple model (|3.3| ), ( |3.4D would 
have to be augmented by terms describing an interaction of instantons and 
anti-instantons M, Esll. 
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Figure 6: /9i -dependence of (Q^)- The numerical data are represented by 
filled circles and we connect the symbols to guide the eye. The numbers on 
the top of the figure give our results for the lattice spacing. 



Let us now analyze the /3i-dependence of (Q^). It is expected that the 
chiral susceptibility x = {Q"^) drops considerably as one goes from the 
confining, chirally broken to the deconfined, chirally symmetric phase. In 
Fig. I we show our results for (Q^) at /?i = 8.10, 8.20, 8.30 and 8.45. Ah data 
are from 6 x 16^ lattices such that after a simple rescaling of the vertical 
axis. Fig. ^ is also a plot for the topological susceptibility. It is obvious that 
with increasing (3i the values of (Q^), and thus the values of x, drop consid- 
erably. This implies that isolated topological excitations are much rarer in 
the high temperature phase = 8.20, 8.30, 8.45) and that their abundance 
decreases further as one moves away from the critical jSi . This fits nicely into 
the above discussed picture where the interaction of topological excitations 
gives rise to the non-vanishing eigenvalue density at the origin p{0) in the 
low temperature phase. Conversely, in the high temperature phase there 
are simply not enough topological excitations to create a sufficient number 
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of interacting pairs which then would build up a non-vanishing eigenvalue 
density /j(0). Of course, the observable (Q^) is a measure for the abundance 
of zero modes, while a non-vanishing p(0) is built up from the near zero 
modes. However, it is quite plausible to assume that as the net number 
of isolated topological excitations decreases also the number of interacting 
pairs of topological excitations goes down. In the next section we will con- 
firm this expectation by analyzing an observable which is directly sensitive 
to the number of interacting pairs. 

4 Results for the eigenvectors 

4.1 A qualitative discussion of some examples 

Before we show our results for the eigenvectors let us first introduce some 
notation. Eigenvectors will be denoted by V- They obey the eigenvalue 
equation 

Dip = Xip. (4.1) 

The eigenvectors have space, color and Dirac indices which we denote by x, c 
and d respectively. For each combination of these indices, the corresponding 
entry il!{x, c, d) is a complex number. We define a scalar density p{x) and a 
pseudoscalar density p5(x) by 

P(^) = Y '>P{x,c,d)* 'tp{x,c,d) 

c,d 

P5{x) = ^ix^c,d)* {^5)^4' 'tp{x,c,d') . (4.2) 

c,d,d' 

Here * denotes complex conjugation. Note that the summation over the 
color index makes the densities in Eq. ( |4.2| ) gauge invariant. Due to the 
normalization tp^tp = 1 of the eigenvectors the scalar density sums up to 1, 
i.e. 

Yp{x) = 1. (4.3) 

X 

In Fig. 1^ we show x-y time-slices for the pseudoscalar densities for two 
characteristic situations (both eigenvectors are for configurations from the 
/?! = 8.45 ensemble on 6 x 16^ lattices). A similar behavior for the eigen- 
vectors of the staggered Dirac operator was observed in [|14|. A study of 
staggered Dirac operator eigenmodes in finite temperature SU(2), |]l^], also 
found localized modes. 
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Figure 7: The 6 x-y time-slices of the pseudoscalar density P5{x) for a zero 
mode (top left plot) and a configuration with an interacting caloron anti- 
caloron pair, a so called topological molecule (bottom right plot). The data 
are from the /3i = 8.45 sample on the 6 x 16^ lattice. 
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In the top left plot we display the six (the temporal extent of our lat- 
tice is Lt = 6) x-y time-slices of P5(x) for an eigenvector with a single 
real eigenvalue. As already discussed above, these modes correspond to 
isolated topological excitations which give rise to a zero mode (in our set- 
ting an eigenvalue with exactly vanishing imaginary part). Let us discuss 
some characteristic features of the pseudoscalar density and compare them 
to the properties known for a classical caloron. For a caloron configuration 
the corresponding eigenvector (i.e. the zero mode) roughly traces the action 
density of the gauge configuration, i.e. p{x) is large where F^^^{x)F^^^{x) is 
large. Thus p{x) is localized in space as is the underlying caloron config- 
uration while it is typically more extended in the time direction, where it 
can interact with itself around the short, compactified time extension of the 
finite temperature lattice. In addition the eigenvector for a classical caloron 
is chiral such that p^{x) = ±p(x). The density p^{x) displayed in the top 
left plot of Fig. shows exactly these features of a zero mode corresponding 
to a caloron field. The density P5{x) is localized and positive, i.e. in this 
particular example it corresponds to a caloron which has positive chirality 
(negative chirality, i.e. P5{x) < would correspond to an anti-caloron) . 

In configurations with several calorons we see several non-degenerate 
real eigenvalues, each having an eigenvector similar to that of the top left 
plot of Fig. ^. Because the modes are not degenerate, we do not observe 
strong mixing of the eigenvectors, but instead each of them is essentially 
tracing only one of the topological excitations. For a Dirac operator which 
is an exact solution of the Ginsparg- Wilson equation one has exact zero 
eigenvalues. So in e.g. a two caloron configuration the two zero modes would 
mix and p{x) and ^5(2;) would in general show peaks at the positions of both 
of the calorons. 

Let us now turn to an eigenvector which has an eigenvalue with a rela- 
tively small but non-zero imaginary part. To be more precise, in the example 
shown here from the /?i = 8.45 ensemble this is a mode at the edge of the 
spectral gap. According to the caloron picture these modes are due to a 
pair of a caloron and an anti-caloron which approach each other and start 
to interact, such that the Dirac operator no longer can distinguish them as 
two isolated topological excitations and the corresponding two zero modes 
are split into a complex conjugate pair. We will refer to these configurations 
of a caloron interacting with an anti-caloron as topological molecule. The 
dipole structure of such a topological molecule is nicely illustrated in the 
bottom left plot. The characteristic feature is a positive peak of P5{x) (on 
the t = 3 slice) next to a negative peak (on the t = 6 slice). 

Finally we remark, that when inspecting P5{x) for an eigenvector with 
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an eigenvalue in the bulk of the distribution, i.e. an eigenvalue with a large 
imaginary part one finds that any trace of a localized topological object has 
been washed out and the pseudoscalar density is dominated by quantum 
fluctuations. 



4.2 Localization properties of the eigenvectors 

In the last section we have presented characteristic examples of eigenvec- 
tors which display some of the features that are expected from models of 
topological excitations. In order to go beyond illustrating these properties 
merely by examples, in this section we now analyze the sample behavior of 
localization properties of the eigenvectors. 

A convenient observable for the localization of an eigenvector is the so- 
called inverse participation ratio which is widely used in solid state physics. 
It is defined by 

I = VJ2 Pi^f ■ (4-4) 

X 

Here V denotes the four- volume LtL^. From its definition in Eq. ( [4.2| ) 
it follows that p{x) > for all x. Taking into account the normalization 
Eq. ( [4.3| ) one finds that a maximally localized eigenvector which has support 
on only a single site x' must have p{x) = Sx^x'- Inserting this into Definition 



(4.4) for the inverse participation ratio one finds that a maximally localized 
eigenvector has I = V. Conversely, a maximally spread eigenvector has 
p{x) = 1/V for all X. In this case, the inverse participation ratio gives a value 



of / = 1. Similarly to the scalar inverse participation ratio I of Eq. (4^) we 
can also define a measure for the localization of the pseudoscalar density 
P5{x) (see Eq. (P)), 

h = VJ2p5{x)^. (4.5) 

X 

We will refer to Is as the pseudoscalar inverse participation ratio (compare 



[13 1 for an alternative observable sensitive to the localization of the pseu- 
doscalar density). 

In Fig. ^ we show the probability distribution P{I) of the inverse par- 
ticipation ratio. We present the results for all lattice sizes of the /?i = 8.10 
and Pi = 8.45 samples, where for the latter case we again distinguish be- 
tween configurations with a real Polyakov loop and configurations in the 
complex sector. In order to avoid the problem of different cut-off values for 
the spectrum when keeping the number of computed eigenvalues fixed at 50 
for all lattice sizes (compare the discussion in Section 2.3), the samples for 
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Figure 8: Distribution of the inverse participation ratio. 



different volumes were cut off at a common value of the imaginary part of 
the eigenvalues. 

A general feature for all plots in Fig. ^ is the increasing probability of 
finding localized objects (large values of /) as the volume gets larger. This 
is consistent with the property of models of topological excitations which 
also grow in number as the volume increases. When comparing the distri- 
bution in the low temperature phase = 8.10) to the sector of the high 
temperature phase = 8.45) with real Polyakov loop we do not see a large 
effect. This means, that the molecules of calorons which form in the chirally 
symmetric phase still are relatively localized objects. More significant is the 
change when one compares the plots for the real and complex sector of the 
Polyakov loop in the chirally symmetric phase {Pi = 8.45, the two plots at 
the bottom of Fig. |8|). In the complex sector we find the modes to be much 
less localized, i.e. the distribution P{I) is much narrower and its peak is 
shifted to smaller values of /. An equivalent dependence of the inverse par- 
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Figure 9: The inverse participation ratio as a function of |Im A|. 



ticipation ratio on the Z3 sector was observed in a study with the staggered 
Dirac operator in |14]. The dependence of / on the phase of the Polyakov 
loop can be understood from a property of the classical solutions. Also for a 
classical caloron one finds that the eigenvector becomes more spread out as 
one performs a Z3 transformation of the gauge field which rotates the phase 
of the Polyakov loop (see Appendix B). 

It is also quite instructive to plot the inverse participation ratio of an 
eigenvector against its eigenvalue. In Fig. ^ we show such a plot. We binned 
the imaginary parts Im A of the eigenvalue spectrum and for each bin com- 
puted the average size (/) of the inverse participation ratio for the corre- 
sponding eigenvectors. Eigenvectors which have a real eigenvalue, i.e. the 
zero modes, were left out. Let us first discuss the plot for the chirally broken 
phase = 8.10, top plot in Fig. ^). It is obvious that for all lattice sizes, 
the most localized states (they have the largest values of /) sit near Im A = 
0. This behavior can again be understood in the context of the picture of a 
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fluid of topological excitations. As already discussed in the last section, the 
near zero modes are believed to come from pairs of topological excitations 
which overlap slightly, i.e. large topological molecules. As long as the two 
excitations are relatively remote, the corresponding eigenvector is still rela- 
tively localized, i.e. its I is large. The closer the two partners come the more 
of their localization is washed out and the corresponding eigenvalue moves 
away from the origin along the imaginary axis. This behavior is exactly 
reflected in the plot where the most localized states are near the origin and 
the localization becomes washed out as one moves along the imaginary axis. 

When turning to the chirally symmetric phase {f3i = 8.45) the behavior 
is quite similar. The only difference is the appearance of the gap such that 
for values of |Im A| below the edge of the gap we find no modes at all. Above 
the edge, however, we again observe that the most localized states are near 
the edge and the localization becomes washed out as one moves further away 
up or down the imaginary axis. For the sector with real Polyakov loop (left 
hand side plot in the bottom row of Fig. ^) this behavior is quite pronounced, 
while due to the spreading effect of a Z3 rotation discussed already above 
and in Appendix B, the signal is much weaker in the sector with complex 
Polyakov loop (right hand side plot in the bottom row). 

In Section 3 we have seen that (Q^) drops considerably as one increases 
/?! and goes over to the high temperature phase. (Q^) is a measure for the 
net number of isolated topological excitations. In order to study topological 
molecules, which do not give rise to zero modes, we analyze the probabil- 
ity for encountering localized modes. We define N{c) to be the number of 
eigenvectors which are not zero modes (the corresponding eigenvalue is com- 
plex and ip^^5ip vanishes) and have an inverse participation ratio larger than 
some cut c. Equivalently we define A^5(c) to be the number of eigenvectors 
with a pseudoscalar inverse participation ratio Is > c. 

In Fig. we show our results for the averages {N{c)) and {N^{c)) for 
three different values of c. All data were computed on lattices of size 6x16^ 
and the values of f3i are /3i = 8.10,8.20,8.30,8.45. 

It is obvious, that for all three values of c there is a clear drop of the 
averages of N(c) and N^{c) for the numbers of eigenvectors with I > c and 
I5 > c respectively, as one increases Pi. As noted above, zero modes were 
omitted in the evaluation of iV(c) and N^ic) such that the observable does 
not include isolated zero modes and is sensitive only to pairs of topological 
excitations which, according to the picture of topological molecules build 
up the eigenvalue density p{0) at the origin. The drop of N{c) and N^{c) 
indicates that the abundance of the topological molecules goes down as 
Pi increases. In the high temperature phase their number has become so 
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Figure 10: Average of N(c) and A^5(c) the numbers of eigenvector with I > c 
and I5 > c respectively, as a function of Zero modes are omitted in the 
evaluation of N{c) and Nq{c). We display our results from 6 x 16^ lattices 
for three values of the cut c. 



small that they no longer create a nonvanishing value of p{0) and the chiral 
condensate vanishes. 

It is interesting to note that N{c) and Nr^{c) drop at essentially the same 
rate. Furthermore this holds individually for both the real and the complex 
sector of the Polyakov loop even though the shape of the drop is slightly 
different in the two sectors. This property of a simultaneous drop of the 
average N{c) and N^{c) further supports the above interpretation that the 
topological molecules simply become more dilute with increasing Pi (which 
causes the drop of /) but keep their local chirality (which leads to a decrease 
of N5 at the same rate as the decrease of N) . 



25 



4.3 Local chirality of the eigenmodes 



In a recent publication [|16[ Horvath, Isgur, McCune and Thacker proposed 
testing the instanton anti-instanton picture of zero temperature QCD by 
analyzing the local chirality of the eigenmodes of the Dirac operator. If the 
near zero modes correspond to interacting instantons and anti-instantons, 
the corresponding peaks in p{x) (compare Eq. ( |4.2| )) should be either left 
or right handed. An example of this behavior is shown in Fig. where left 
handed peaks show up as positive values of p5(x) and right handed peaks are 
negative. However, using the Wilson Dirac operator and a relatively small P 
the authors of |]l6| did not observe local chirality of their near zero modes. In 
the meantime their observable has been analyzed by different groups using 



slightly different settings [17|-|2C] and all those later publications do indeed 
see local chirality supporting the instanton anti-instanton picture. 

To clarify this point we also analyzed the local chirality variable of p^ . 
Before we present our numerical results we briefly repeat the definition of 
the local chirality variable of Similar to the densities p{x) and p^ix) of 
Eq. ( |4.2| ) we can now define local densities p+{x) and p_(x) with positive 
and negative chirality 

P±{x) = E Hx,c,dy {P±)d,d' ^{x,c,d') , (4.6) 

c,d,d' 

where P± are the projectors onto positive and negative chirality P± = (1 =b 
75)/2. It is now interesting to analyze locally for each lattice point x the ratio 
pj^{x) /p-.{x). For a classical instanton in the continuum, the corresponding 
zero mode ip{x) has positive chirality and the density p+{x) is positive for all 
X while p- [x) vanishes everywhere. Thus the ratio of the two always gives 
oo. For an anti-instanton the roles of p+ and p- are exchanged and the 
ratio is always 0. When now analyzing the eigenvectors for an interacting 
instanton anti-instanton pair this property should still hold approximately 
near the peaks. The ratio p+{x) /p-{x) is expected to be large for all x near 
the instanton peak of the pair and small for all x near the anti-instanton 
peak. In a final step Horvath et al. map the two extreme values oo and of 
the ratio pj^{x) /p-.{x) onto the two values —1 using the inverse tangent. 
One ends up with the local chirality variable X{x) defined as 



When the eigenvectors for the near zero eigenvalues correspond to instanton 
anti-instanton pairs then the values of X{x) should cluster near +1 and —1 
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when one chooses lattice points x near the peaks of the scalar density p{x). 
One can choose different values for the percentage of points x shown and 
here we will present results for cuts of 1%, 6.25% and 12.5%. This means 
that we average over those 1% (6.25%, 12.5%) of all lattice points which 
support the highest peaks of p{x). The smallest cut-off of 1% will give the 
best signature since only the highest peaks which are not so much affected 
by quantum fluctuations do contribute. On the other hand such a small 
cut-off may not be very conclusive, since typically instanton models have a 
packing fraction of instantons considerably larger than 1% (see ||^, ^]). 

We would like to remark that in |18| a modification of the local chirality 
variable X was presented, which uses a bi-orthogonal system to construct 
p±{x). This considerably improves the signal for the Wilson Dirac operator. 
We tested the version proposed in [18| and found only a minor improvement 
of our results. This is due to the fact that we use an approximation of a 
Ginsparg- Wilson Dirac operator. A solution of the Ginsparg- Wilson equa- 
tion is a normal operator which can be diagonalized with a unitary transfor- 
mation and the bi-orthogonal system then reduces to the right eigenvectors 
used here. 

We start the discussion of the local chirality with the zero modes. In 
Fig. ^ we show our results for the distribution of X on both sides of the 
phase transition. The data for the chirally broken phase (/3i = 8.10) was 
computed from 55 gauge configurations. For the chirally symmetric phase 
(/?i = 8.45) we also analyzed 55 gauge configurations but here we distinguish 
between configurations with real and complex Polyakov loop, such that the 
sample size for each of the sectors individually is only one third, respectively 
two thirds of the overall size of 55. In addition, zero modes are much rarer 
for the /3i = 8.45 sample (compare e.g. Fig. Q) and the two effects combined 
account for the larger error bars in the (5i = 8.45 data. We also attribute 
the slight asymmetry observed in the data for the complex sector of the 
Polyakov loop to a fluctuation in the relatively small statistics for this plot 
which seems to contain more left-handed than right-handed modes (compare 
also Fig. ^). All data we show were computed for the largest volume 6 x 20^. 
We remark that our data for P{X) are normalized such that the area under 
the curve is 1. 

All three plots in Fig. |ll| show a clear dip of P{X) near X = 0. Further- 
more the shape of P{X) is almost independent of the percentage of lattice 
points chosen. This shows that zero modes of the Dirac operator are locally 
chiral up to at least one eighth of the total volume. It is remarkable that 
the local chirality signal for the zero modes in the chirally broken phase 
(/?! = 8.10) is less pronounced (the two peaks are smaller and wider) than 
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Figure 11: Local chirality for the zero modes. We present data for both 
sides of the phase transition. For the chirally symmetric phase, we distin- 
guish between the real and the complex sector of the Polyakov loop. We 
use different values for the cut-off on the percentage of supporting lattice 
points: 1%, represented by filled circles, 6.25%, open squares and 12.5%, 
filled triangles. The data were computed on the 6 x 20^ lattice. 



28 



for the zero modes in the chirally symmetric phase. This might be due to 
the fact that the smaller value of /3i allows for fluctuations on smaller scales 
(measured by the number of lattice points) which spoil the signal for local 
chirality. These UV fluctuations become more suppressed as one increases 
Pi- Still we can summarize: We find a clear signal for local chirality of the 
zero modes on both sides of the phase transition. 

Let us now turn to the non zero modes in the chirally broken phase 
(/?! = 8.10). In Fig. 12 we show our results. Again they were computed 
from 55 gauge configurations on the 6 x 20^ lattice and we use the same cut- 
off values as for the zero modes, i.e. 1%, 6.25% and 12.5%. In addition we 
binned the eigenvectors used for the computation of the local chirality with 
respect to the imaginary parts of the corresponding eigenvectors. We divided 
all eigenvectors into three bins such that their eigenvalues obey |Im A| < ci 
for the first bin, ci < |Im A| < C2 for the second bin and C2 < |Im A| for the 
third bin. The thresholds q were set such that on the average the first bin 
contained 4 eigenvectors and the other two bins each 8 eigenvectors (for the 
computation of the local chirality we always stored a total of 20 eigenvectors 
which are not zero modes). The values of the q are given by ci = 0.018 
and C2 = 0.057. The first bin (top plot in Fig. p^ ) contains the very near 
zero modes, while in the second and third bin (middle and bottom plots in 
Fig. 12) we give the results for eigenvectors with larger |Im A|. 



For the very near zero modes (|Im A| < ci, top plot in Fig. 12) we find 
a very pronounced signal for the local chirality. In fact it is almost indistin- 
guishable from the result for the zero modes shown in the top plot of Fig. pT| 
(note the different scale on the vertical axis in Fig. 11 and Fig. 12). This 
shows that the topological molecules responsible for the near zero modes are 
still very much reminiscent of the localized structure of the two partners, 
i.e. the caloron and the anti-caloron. Furthermore the dependence for the 
cut-off on the number of supporting points is very weak, i.e. the 1% and 
the 12.5% curves are still very similar. This indicates that these modes are 
relatively spread out, such that, similar to the zero modes, more than one 
eighth of the volume is occupied by excitations with a definite chirality. 

As the eigenvalues move away from the origin, the local chirality of the 
corresponding eigenvectors becomes less pronounced. This can be seen from 
the middle and bottom plots in Fig. ^ where we show our results for the 
second and third bin of |Im A|. Firstly we find that the peaks are less 
pronounced and secondly, that now the curve is much more sensitive to the 
cut-off for the percentage of supporting lattice points. The interpretation 
of both these observations is straightforward. According to the instanton 
picture, eigenvalues higher up in the spectrum are produced by pairs of 
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X 



Figure 12: Local chirality for non zero modes in the chirally broken phase. 
We bin the eigenvectors with respect to the imaginary parts of the corre- 
sponding eigenvalues such that the first bin (|Im A| < ci, top plot) shows 
the very near zero modes, while the other two bins (ci < |Im A[ < C2, middle 
plot and C2 < |Im A|, bottom plot) show the local chirality for eigenvectors 
with eigenvalues higher up in the spectrum. The choice of the thresholds Cj 
is discussed in the text. 
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<L> real <L> complex 




Figure 13: Local chirality for non zero modes in the chirally symmetric 
phase (/3i = 8.45,6 x 20^). We divide our data with respect to the sector 
of the Polyakov loop. The left column of plots contains the results for real 
(L), the right column shows data for complex (L). In addition we bin the 
data with respect to |Im A| such that the first bin (top plots) shows the local 
chirality for eigenvectors with eigenvalues close to the edge of the spectral 
gap, while the other bins (middle and bottom plots) show eigenvectors with 
eigenvalues more remote from the edge. 
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topological excitations which are overlapping more than the excitations for 
the near zero modes. Thus the localization for the two partners in a molecule 
is more washed out and so is the local chirality. This washing out is also 
reflected in the stronger dependence on the percentage for the supporting 
lattice points. If the overall height of the peaks in p{x) is going down, then 
the 1% cut may still see only the highest peaks which are still locally chiral, 
but the larger values 6.25% and 12.5% start to cut also into the quantum 
fluctuations which do not have local chirality. To summarize the discussion 
of Fig. ^ we find good agreement with the picture of topological molecules 
which produce the nonvanishing eigenvalue density near the origin and lead 
to eigenvalues higher up (down) in the spectrum as the two partners in the 
molecule start to overlap more. 

Finally we discuss our results for the non zero modes in the chirally 
symmetric phase (/3i = 8.45,6 x 20^). Fig. ^ shows our data which were 
again divided into the two sectors for the Polyakov loop and the left column 
of plots shows the sector with (L) real and the right column shows the 
complex sector. As in the chirally broken phase we bin the eigenvectors with 
respect to the imaginary parts of the corresponding eigenvalues. The first 
bin (|Im A| < ci, top plot in Fig. 13) contains the data for eigenvectors with 
eigenvalues very close to the edge of the gap, while the other two bins (ci < 
|Im A I < C2, middle plot and C2 < |Im A|, bottom plot) show the results for 
eigenvectors with eigenvalues more remote from the edge. Again we chose 
the values of the thresholds such that the first bin contains an average 
of 4 eigenvectors while the other two bins each contain 8 eigenvectors. This 
binning is done for each of the sectors individually. The resulting values for 
the thresholds at /3i = 8.45 are then given by ci = 0.207 and C2 = 0.230 for 
the real sector and ci = 0.081 and C2 = 0.112 for the complex sector. 

It is obvious, that for the real sector (left column of plots) the signal 
for local chirality has vanished entirely. For all three bins and for all three 
percentages of supporting lattice points we do not see a dip in P{X). This 
shows that eigenvectors with local chirality do not play a substantial role in 
the real sector of the chirally symmetric phase at Pi = 8.45. Similarly in the 
data for the complex sector of the Polyakov loop only the 1% cut-off, i.e. only 
the highest peaks of the scalar density, produces a small dip in P(X), which, 
however, vanishes as one increases the percentage of supporting points to 
6.25% or even 12.5%. Thus only about 1% of the fluctuations in the complex 
sector of the /3i = 8.45 data show local chirality. Based on Fig. IC and the 
subsequent discussion we interpret the failure to find non-zero modes with 
local chirality as further evidence for the diluteness of topological molecules 
in the chirally symmetric phase. 
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5 Summary 



In this article we have presented a large scale study of eigenvectors and 
eigenvalues of the lattice Dirac operator. We analyzed the Dirac operator 
for background gauge fields from quenched finite temperature ensembles on 
both sides of the QCD phase transition. Using three different lattice sizes 
we analyzed the volume dependence of our observables. We used a Dirac 
operator which is an approximate solution of the Ginsparg- Wilson equation 
and has good chiral properties. 

Our observables focused on testing the topological nature of the gauge 
field excitations as seen by the Dirac operator. In particular we analyzed the 
density of the eigenvalues, the topological charge defined through the index 
theorem, scalar and pseudoscalar densities for the eigenvectors, their local- 
ization properties and a recently proposed observable for the local chirality 
of the eigenvectors. Our main results are: 

• For the low temperature phase we find a density p{X) of eigenvalues 
which extends all the way to the origin which, according to the Banks 
Casher relation, gives rise to a non-vanishing chiral condensate in this 
phase. In the high temperature phase we observe a gap in the spectral 
density and thus a restoration of chiral symmetry. 

• The spectral gap opens up for all phases of the Polyakov loop, i.e. when 

we divide our ensemble into configurations with real Polyakov loop and 
configurations where the Polyakov loop has a phase of ±2tt/3 we find a 
gap for both subsets. This shows that the chiral symmetry is restored 
in both these sectors. 

• When analyzing the topological charge Q as defined by the index theo- 
rem we find a probability distribution and values of (Q^) which can be 
described well by a simple model of topological excitations. The drop 
of (Q^) as one increases (5i reflects a thinning out of the topological 
excitations which prevents the formation of a chiral condensate in the 
high temperature phase. 

• Plots of the pseudoscalar densities for a zero mode and a near zero 
mode illustrate the characteristic behavior of the different types of 
excitations. 

• A measure for the localization of the eigenvectors is given by the so 

called inverse participation ratio. For the low temperature phase, we 
found that the most localized eigenvectors have eigenvalues near the 
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origin and the localization of eigenvectors gets washed out as the imagi- 
nary part of the corresponding eigenvalue increases. Similarly one finds 
for the high temperature phase that the most localized eigenvectors 
have eigenvalues near the edge of the distribution. These observa- 
tions support the picture that near zero modes, respectively the near 
edge modes, are the still relatively localized remnants of interacting 
calorons and anti-calorons. 

• In the high temperature phase we find that the eigenmodes in the 
Z3 sector where the Polyakov loop is real are more strongly localized 
than modes in the complex Z3 sector. This behavior supports the 
idea that calorons are relevant, because the classical solution of the 
Dirac equation in a caloron gauge field shows the same behavior, see 
Appendix B. 

• When plotting the number of eigenvectors with a large localization as a 
function of Pi we found a drop when increasing /3i , similar to the drop 
observed for (Q^). This shows that also interacting pairs of calorons 
and anti-calorons become less abundant at larger Pi and in the high 
temperature phase their number is no longer sufficient to build up a 
non-vanishing chiral condensate. 

• Finally we analyzed the local chirality of the eigenvectors at different 
values of the cut-off for the scalar density. For the chirally broken 
phase our results for the near zero modes confirm the pattern of local 
chirality as expected for molecules of topological excitations. In the 
chirally symmetric phase we still see modes with local chirality near 
the edge of the gap, but they are much more diluted by bulk modes 
dominated by quantum fluctuations. 

Our results strongly support a picture of interacting topological excita- 
tions which are responsible for the breaking of the chiral symmetry at low 
temperatures, while in the high temperature phase they are no longer suffi- 
ciently abundant and the chiral symmetry is restored. 
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A Detailed specification of the chirally improved 
Dirac operator 

In this appendix we describe in more detail the terms in our Dirac operator 
and give the values for the coefficients which we use for the four ensembles 
of quenched gauge field configurations. 



Clifford generator Generating path Name of coefficient 



I 


<> 


Sl 


I 


< 1 > 


S2 


I 


< 1,2 > 


S3 


I 


< 1,2,3 > 


S5 


I 


< 1,1,2 > 


S6 


I 


< 1,2,-1 > 


S8 


I 


< 1,2,3,4 > 


SlO 


I 


< 1,2,-1,3 > 


fill 


I 


< 1,2,-1,-2 > 


■S13 


71 


< 1 > 


Vl 


71 


< 1,2 > 


V2 


71 


< 1,2,3 > 


Vi 


71 


< 2,1,3 > 


V5 


71 72 


< 1,2 > 


tl 


71 72 


< 1,2,3 > 


t2 


71 72 


< 1,3,2 > 


t3 


71 72 


< 1,2,-1 > 


t5 


71 72 


< 1,2,-1,-2 > 


tl5 


75 


< 1,2,3,4 > 


Pi 



Table 3: Description of the terms in our D. 

As has been pointed out in Section 2.1, the most general Dirac operator 
D can be expanded in the series (2^). Each term in this series is char- 
acterized by three pieces: A generator of the Clifford algebra, a group of 
paths and a real coefficient. The paths within a group can have different 
signs which are determined by the symmetries, C, P, 75-hermiticity and ro- 
tation invariance. The symmetries also determine which paths are grouped 
together. Thus it is sufficient to characterize a group of paths by a single 
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generating path and all the other paths in the group as well as their relative 
sign factors can be determined by applying the symmetries. 

In addition, for the vector and tensor terms appearing in our D it is 
sufficient to give the paths only for one vector (tensor) since rotation invari- 
ance immediately fixes the structure for the other vector (tensor) terms. In 
order to describe our D we start with listing the three determining pieces 
for each term in Table ^. 

In Table ^ we list the values of the coefficients Si,Vi, U and pi for the 
different values of the inverse gauge coupling (5i for the ensembles used in 
this work. 





Pi = 8.10 


Pi = 8.20 


Pi = 8.30 


Pi = 8.45 


Sl 


+1.54498 


+1.54590 


+1.54688 


+1.54737 


S2 


-0.06169 


-0.06063 


-0.05997 


-0.05892 


S3 


-0.01448 


-0.01449 


-0.01448 


-0.01451 


Sb 


-0.00262 


-0.00258 


-0.00255 


-0.00251 


S6 


+0.00220 


+0.00215 


+0.00210 


+0.00206 


ss 


-0.00540 


-0.00534 


-0.00532 


-0.00525 


SlO 


-0.00053 


-0.00052 


-0.00051 


-0.00050 


Sll 


-0.00118 


-0.00117 


-0.00117 


-0.00116 


Sl3 


+0.00780 


+0.00778 


+0.00777 


+0.00775 


Vl 


+0.10975 


+0.11083 


+0.11063 


+0.11240 


V2 


+0.01770 


+0.01724 


+0.01707 


+0.01654 




+0.00744 


+0.00767 


+0.00784 


+0.00804 




+0.00182 


+0.00193 


+0.00198 


+0.00209 


h 


-0.09874 


-0.09862 


-0.09860 


-0.09843 


t2 


-0.00309 


-0.00303 


-0.00299 


-0.00292 


t3 


+0.00227 


+0.00224 


+0.00221 


+0.00218 


t5 


-0.00663 


-0.00651 


-0.00641 


-0.00627 


tl5 


-0.00361 


-0.00358 


-0.00356 


-0.00354 


Pi 


-0.00911 


-0.00907 


-0.00905 


-0.00901 



Table 4: The numerical values of the coefficients Si,Vi, ti and pi for the 
different ensembles of gauge fields. 
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B Quark zero modes 



B.l Solutions of the field equation 

In this section we discuss the quark zero modes in the caloron gauge field so- 
lutions of Harrington and Shepard [46|. In particular, we want to investigate 
their dependence on the Z3 sector. 

't Hooft has shown that one way of generating multi-instanton solutions 
of the gluon equations of motion is to consider gauge fields of the form 

A^^^,[x) = ia'Z,d,\np{x), (B.l) 

where c and c' are color indices and the a are a set of 2 x 2 color matrices. 
A suitable basis for the a is given in |47]. After a little algebra one finds 
that this A field satisfies the equations of motion if 

d'^p{x) = (B.2) 
or d'^p{x) oc p^{x) . (B.3) 

We know from the Atiyah-Singer index theorem that gauge field configura- 
tions with topological charge must have fermion zero modes. Following [^] 
these modes can be written in a form analogous to (|B.1|), 



^^,{x) = p'/\x)a^^,d,(^^^ , (B.4) 
where c is a color index, and d a Dirac index. The relation of the a matrices 



to the a matrices, and a choice of a is given in |47|. The ansatz ( |B.4| ) gives 
a zero mode of the Dirac operator if 

p{x)d\{x)-x{x)d^pix) = 0. (B.5) 

We now need a solution of Eq. ( |B.2| ) that also satisfies the periodic 
boundary condition in the time direction, 

pir,t + l/T)=pir,t) , (B.6) 

where T is the temperature (we have set Boltzmann's constant to 1). A 
simple solution of this type, IM] , is 



V TrTi?^ sinh(27rTr) 

^k^^ r2 + (t - k/T)^ ~ ^ r [cosh(2^rr) - cos(2^rt)] ^ ^ 
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where = rf + + r^. This sokition corresponds to a caloron of radius R 
sitting at the origin. We can use translation invariance to give solutions at 
other locations, and add together solutions of the above form to get multi- 
caloron solutions. 

We now look for a fermion zero mode in the solution ( |B.7| ) , which should 
obey anti-periodic boundary conditions 

^{r,t + l/T) = -^{r,t) . (B.8) 

■0 depends linearly on x, so x must obey the same boundary conditions as 
ip. The function 



1 ^ R 



TR cos(7rTt) sinh(7rrr) ^ 



r [cosh(27rrr) — cos(27rTt)] 



satisfies eq. ( |B.5f ), and has the correct behavior under i — > t + 1/T. This 
has the same form as the general multi-instanton solution given in 

The solution ( [B.9D corresponds to the real sector of the Z3 symmetry. We 
are also interested in knowing what the zero-mode looks like in the complex 
sectors. There are two possible approaches to this problem. 



The direct approach is to use the fact that eq. (BT) gives an SU(2) A- 
field. When we embed this solution in SU(3) the ^-field has components 
proportional to the three generators and of SU(3). The generator 

commutes with these three matrices, so we still have a solution of the 
equation of motion if we add an arbitrary constant A field in the 8-direction 
of SU(3). This can be used to produce a caloron solution with an arbitrary 
Polyakov loop value. See fi^ for methods of constructing SU(2) solutions 
with non-trivial Polyakov loop. 

A simpler method is to note that a Polyakov loop can be gauge trans- 
formed to look just like a fermion boundary condition of the form 

i^{f,t + l/T) = -e'\{r,t) . (B.IO) 

For the complex Z3 sectors of high-temperature QCD we are interested in 



9 = ib27r/3. A solution of eq. (|B.5| ) that obeys Eq. ( B.IO ) is given by 



1 °° R 



e TR [cosh(0Tr) cos{TTTt) sinh(7rrr) — i sinh(0Tr) sin(7rTt) cosh(7rTr)] 

r [cosh(27rTr) — cos(27rrt)] 
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Figure 14: The scalar density p{x) for a quark zero mode in the real Z3 
sector (upper plot) and the complex Z3 sector (lower plot). In both cases 
the caloron radius is 1/T. 
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valid for —it < < it. At large r the probability density found from 
eq. ( B.ll| ) drops off like r^^exp[— 2(7r — |0|)Tr]. This is much faster in 



the real sector (6 = 0) than in the complex sectors {9 = ±2tt/3). 
B.2 Phenomenological implications 

To illustrate the difference between the Z3 sectors, we plot the density of 
the fermion zero mode in a caloron of radius R = 1/T, see Fig. |l^. We show 
a two-dimensional x-t slice through the center of the mode. Both modes 
have been normalized so that the 4-dimensional integral of p is 1. One sees 
clearly that the mode in the real sector is more strongly localized. 



RT 



Figure 15: The inverse participation ratio I of the quark zero- mode cal- 
culated for a caloron in the real Z3 (solid line) and in the complex sector 
(dashed line). For comparison the result for an instanton is shown as the 
dotted line. R is the radius of the caloron and T the temperature. 



A systematic way to compare the localization of the modes (|B.9|) and 



( B.ll ) is to define an inverse participation ratio analogous to that of Section 
4. The scalar density p{x) is defined in exactly the same way in the contin- 
uum as on the lattice, Eq. (^.2]). In the definition of / we simply replace the 



sum in Eq. (^^) by an integral: 

I/V = I (fx p{xf (B.12) 
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where the eigenmode has been normaUzed so that / d'^x p{x) = 1. 

In Fig. 15 we compare the locahzation of the quark zero modes as the 
caloron radius R is varied. When RT <C 1 thermal effects are small, so a 
caloron looks very like an instanton, and the / of the zero mode has only 
a slight dependence on the Z3 sector. This contrasts with the situation 
when the caloron radius is comparable with the inverse temperature, RT ~ 
1. Here the Z3 sector has a dramatic effect, with localization being much 
stronger in the real sector than in the complex. This is reminiscent of the 
effects we see in the high temperature phase (see Section 4). 
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